d = 5;
start = 0;
thisend = 10;

x1 = [start:.1:thisend];
y1 = f(x1);

plot(x1,y1);

quad("f", start, thisend)